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Abstract. The Q./Q. ratio originating from string decays is predicted to be larger 
than unity in proton-proton interaction at SPS energies. The anti-omega dominance 
increases with decreasing beam energy. This surprising behavior is caused by the 
combinatorics of quark-antiquark production in small and low-mass strings. Since this 
behavior is not found in a statistical description of hadron production in proton-proton 
collisions, it may serve as a key observable to probe the hadronization mechanism in 
such collisions. 
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1. Introduction 

Hadron yields and their ratios stemming from the final state of ultra-relativistic heavy- 
ion collisions have been extensively used to explore the degree of chemical equilibration [ 
|l[ || fj], g, |f| 0, H |To| l and to search for evidence for exotic states and phase transitions 
in such collisions [|l|]. Under the assumption of thermal and chemical equilibrium, fits with 
a statistical (thermal) model have been used to extract bulk properties of hot and dense 
matter, e.g. the temperature and chemical potential at which chemical freeze-out occurs [ 

IB BIB; 

The application of a statistical model to elementary hadron-hadron reactions was first 
proposed by Hagedorn [ MJj in order to describe the exponential shape of the ra r -spectra of 



produced particles in p+p collisions. Recent analyses [ 12 1 on hadron yields in electron- 
positron and proton-proton interactions at several centre-of-mass energies have shown that 
particle abundances as well can be described by a statistical ensemble with maximized 
entropy. In fact, the abundancies are consistent with a model assuming the existence of 
equilibrated fireballs at a temperature T w 160 — 170 MeV. These findings have given re- 
newed rise to the interpretation that hadronization in elementary hadron-hadron collisions 
is a purely statistical process, which is difficult to reconcile with the popular dynamical 
picture that hadron production in pp collisions is due to the decay of color flux tubes [ pj| ]. 
In this article we argue that the Q./Q. = Cl + /Or ratio in elementary proton-proton 
collisions is an unambiguous and sensitive probe to distinguish between particle production 
via the breakup of a color flux tube from statistical hadronization [ [l4| ]. 

2. (Anti-)baryon production in sting models 

Color flux tubes, called strings, connect two SU(3) color charges [ 3 ] and [ 3 ] with a linear 
confining potential. If the excitation energy of the string is high enough it is allowed to 
decay via the Schwinger mechanism [ [T^ ], i.e. the rate of newly produced quarks is given 
by: 

exp -TCmi/K (1) 

dp ± 



where K is the string tension and m± — \J p\+ m 2 is the transverse mass of the produced 
quark with mass m. 

However, specific string models may differ in their philosophy and the types of strings 
that are created: 

• In UrQMDf [l(J| the projectile and target protons become excited objects due to the 
momentum transfer in the interaction. The resulting strings, with at most two strings 
being formed, are of diquark-quarktype. 



In NeXuS[ |l7j], the pp interaction is described in terms of pomeron exchanges or 
ladder diagramms. Both hard and soft interactions happen in parallel. Energy is 
shared equally between all cut pomerons and the remnants. The endpoints of the 



Probing Hadronization with Strangeness 



3 



cut pomerons (i.e. the endpoints of the strings) may be valence quarks, sea quarks, 
antiquarks or gluons. 

• In PYTHIA[ |l8|], a similar scheme as in UrQMD is employed. However, hard inter- 
actions may create additional strings from scattered sea quarks. Most strings are also 
of diquark-quark form. 

Fig.|l]a depicts the anti-baryon to baryon ratio at midrapidity in proton-proton interac- 
tions at 160 GeV. The results of the calculations by NeXuS, UrQMD and PYTHIA, which 
are the best established string-fragmentation models for elementary hadron-hadron inter- 
actions, are included in this figure. In all these models, the B/B ratio increases strongly 
with the strangeness content of the baryon. For strangeness |s| =3 the ratio significantly 
exceeds unity. In UrQMD and PYTHIA the hadronization of the diquark-quark strings 
leads directly to the overpopulation of £2. In NeXuS however, the imbalance of quarks and 
anti-quarks in the initial state leads to the formation of q m \ — s sea strings, (the s m \ — <7 sea 
string is not possible). These strings than result in the overpopulation of Q.'s. 



Model I 
Model II 
NA49 



; i: 

1 2 3 12 3 

strangeness s strangeness s 

Fig. 1. Left: anti-baryon to baryon ratio at midrapidity in pp interactions at 160 GeV as 
given by PYTHIA, NeXuS and UrQMD as well as statistical model. Right: anti-baryon 
to baryon ratio in 471 for the same reaction as given by statistical models. Stars depict 
preliminary NA49 data for the B/B ratio at midrapidity. 




3. (Anti-)baryon production in statistical models 

In Fig.|] the string model results are compared with the predictions of statistical models 
(SM). Within the SM, hadron productions is commonly described using the grand canon- 
ical (GC) partition function, where the charge conservation is controlled by the related 
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chemical potential. In this description a net value of a given U(l) charge is conserved on 
average. However, in the limit of small particle multiplicities, conservation laws must be 
implemented exactly, i.e., the canonical (C) ensemble for conservation laws must be used 
[ [l2|, |l9[ Ej. The conservation of quantum numbers in the canonical approach severely 
reduces thermal phase space available for particle production. Thus, exact charge con- 
servation is of crucial importance in the description of particle yields in proton induced 



processes and in e + e~ [ 12|, as well as in peripheral heavy ion collisions [ [20[ ]. 

In Fig.|l] the predictions of two different canonical statistical models for B/B ratios 
in pp collisions are included. The main difference between these models is contained in 
the implementation of baryon number and isospin conservation as well as how additional 
strangeness suppression is introduced. 

(I) The calculation in this statistical model [ [24j is a full canonical one with fixed 
baryon number, strangeness and electric charge identical to those of initial state. Also, an 
extra strangeness suppression is needed to reproduce the experimental multiplicities. This 
is done by considering the number of newly produced (ss) pairs as an additional charge 
to be found in the final hadrons. The ss pairs fluctuate according to a Poisson distribution 
and its mean number is considered a free parameter to be fitted [ |4j. The parameters 
used for the prediction of £1 + /Qr ratio (T, global volume V sum of single cluster volumes 
and (ss)) have been obtained by a fit to preliminary NA49 pp data [ E5|] yielding T = 
183.7 ±6.7 MeV, VT 3 =6.49 ±1.33 and (ss) =0.405 ±0.026 with a yf/dof = 11.7/9. 
It must be pointed out that £2 + /£2~ ratio is actually independent of the (ss) parameter and 
only depends on T and V (see also Fig.|j). 

(II) Here we first approximate the conservation of baryon number and electric charge 
by the GC ensemble. Under thermal conditions at top SPS energies this approximation 
leads to deviations from the exact C results in pp collisions by at most 20-30% [ \2^\ . 
Strangeness conservation is, however, implemented on the canonical level following the 
procedure proposed in [ |2(]] . It accounts for strong correlations of produced strange parti- 
cles due to constraints imposed by the locality of the conservation laws. In pp collisions 
strangeness is not distributed in the whole volume of the fireball but is strongly correlated. 
A correlation volume parameter Vb = AtzRq/3 is introduced, where Rq ~ 1 fm is a typical 
scale of QCD interactions. Previous analysis of WA97 pA data yields: Rq ~ 1 . 12 fm corre- 
sponding to Vb ~ 5.8 fm 3 . Note that hidden strange particles are not canonically suppressed 
in this approach. Analysis of experimental data in AA collisions has shown that T and /jb 
are almost entirely determined by the collision energy and only dependent weakly on the 
number of participants 4jt results of NA49 on p/% and %/A part ratios in pp and PbPb 
collisions coincides within 20-30%. In terms of the SM this can be understood if T and /jb 
in pp and PbPb collisions have similar values. We take T ~ 158 MeV and /jb — 238 MeV 
as obtained from SM analysis of full phase-space NA49 Pb-Pb data [ |l^j. The volume of 
the fireball V ~ 17 fm 3 and the charge chemical potential in pp was then found to reproduce 
the average charge and baryon number in the initial state. 

The predictions of the statistical models are shown in fig [jjright). In these approaches 
the B/B ratio increases linearly with the strangeness content of the baryon. For comparison, 
both figures include preliminary data on the B/B ratio at midrapidity by NA49 [ ^5[| (stars). 
Roughly 45 £2 have been so far extracted, no £2 has been observed. With a 95% confidence 



Probing Hadronization with Strangeness 



7 
6 

IT) 5 

b 

x 4 

P 3 



, UrQMD 


....... n + 




n . 


« — t 




J 


•\'-. 

\ 



-3 



1 



NeXus 


Q + 




n" . 




s — •, 


/ 


¥ 


.J . 


. V 


-2 -1 


1 2 




y 



PYTHIA 



-1 

y 

Fig. 2. Rapidity density of anti-Omegas and Omegas in pp interactions at 160 
predicted by UrQMD, NeXus, PYTHIA. 
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level, the £2/£2 ratio in this measurement is below 0.5 [ p3[ ]. However, it must be stressed 
that in the statistical models one can calculate particle production only in full phase space 
and all quoted predictions refer to fully integrated particle spectra. 

Large deviations of statistical models from the data seen in Fig. lb are to be expected 
as at midrapidity B/B ratios are known to be much larger than in the whole phase space. 
To make predictions for particle ratios at midrapidity one can use the canonical model II. 
Following the procedure described in [ ^] we first choose T ~ 168 MeV, from the fit to 
the midrapidity WA97 PbPb data [ 0], and /j b — 130 MeV to reproduce p/p ~ 0.22 in pp 
collisions. With these parameters the agreement of the model II with NA49 data is seen in 
Fig.l to be quite satisfactory. 



4. Rapidity-, mass- and volume dependencies 

The rapidity dependence of Q. and Q. yield is studied in Fig. || within different string models. 
The results were calculated in pp interactions at 160 GeV within PYTHIA, NeXuS and 
UrQMD (from top to bottom). As can be seen, the ratio is largest around mid-rapidity. 

The Q./Q. ratio is fairly robust - different string-model implementations (PYTHIA, 
UrQMD, NeXuS) all agree in their predictions within ±20%. The value of H/D. > 1 in 
p+p reactions at the SPS is a generic feature of the string-fragmentation. However, as 
shown in Table I the total 4ft-yields of £2 s and £l's may differ by a factor of four among the 
different string models. The statistical models are in general giving more consistent results, 
however deviations up to 20% are not excluded. 

In string models, the particle abundances depend on the parameters chosen for the 
fragmentation scheme, while in statistical models they reflect the differences between the 
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Fig. 3. Fragmentation of a color field into quarks and hadrons. While in large strings £2s 
and £2s are produced in equal abundance (top), small strings suppress £2s at the string ends. 



ensembles chosen. Thus, the absolute yields allow to distinguish between the implementa- 
tions once experimental data becomes available. 

In order to understand the large Q./Q. values predicted by string models one elucidates 
in Fig.|| the color flux tube break-up mechanism. Fig.|| shows the fragmentation of the 
color field into quark-antiquark pairs, which then coalesce into hadrons. While in large 
strings Q. & and £2s are produced in equal abundance (a), low-mass strings in UrQMD 
suppress £2 production at the string ends (b), while in NeXuS £2s are enhanced (c). Thus, 
the microscopic method of hadronization leads to a strong imbalance in £2/£2 ratio in low- 
mass strings. 

The £2/£2 ratio depends in a strongly non-linear fashion on the mass of the fragmenting 
string. Fig.|] shows the Q./Q. ratio as a function of the mass of the fragmenting string (i.e. 
different beam energies in pp). One clearly observes a strong enhancement of £2 production 
at low energies, while for large string masses the ratio approaches the value of Q./Q. = 1 
(which should be reached in the limit of an infinitely long color flux tube). 

However, it should be noted that recently a new class of string models utilizing parton- 
based Gribov-Regge theory has been proposed which are capable of generating an B/B 
ratio of less than one [ |2(]]. 

Statistical models, on the other hand, are not able to yield a ratio of £2/£2 > 1. This 
can be easily understood in the GC formalism, where B/B ratio is very sensitive to the 
baryon-chemical potential /jr- For finite baryon-densities, the B/B ratio will always be < 1 
and only in the limit of — may £2/£2 = 1 be approached. These features survive in the 
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Model 



n(xi(r 4 ) n(xi(r 4 ) 



NeXus 0.48 0.79 

PYTHIA 0.17 0.30 

UrQMD 0.66 1.05 

Canonical Model I 0.46 0.31 

Canonical Model II 0.41 0.24 



Table 1. 4tc particles yields in pp collisions at 160GeV 
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Fig. 4. Left: Q./D. ratio as a function of string mass. Right: £2/£2 ratio as a function of the 
volume in Modell II. 



canonical framework, where the GC fugacities are replaced by ratios of partition functions 
[ O]. This is shown in Fig.|](right) where the ratio £2/£2 in pp collisions (according to 
the previously described model I) is plotted as a function of volume for four different 
temperatures. Hence, finite size corrections in the statistical model actually lead to the 
opposite behavior[ ^TJ in the ratio of D./D. vs. system-size (i.e. volume replacing string- 
mass) than that observed in the fragmenting color flux tube picture. 

5. Balance Functions 

While the £2 /Q ratio provides us with important information on the dynamics of hadroniza- 
tion, it does not yield any information on the time-scales at which hadronization occurs: 
Balance functions offer a unique model-independent formalism to probe the time-scales of 
a deconfined phase and subsequent hadronization [ ^7j|. Late-stage production of quarks 
could be attributed to three mechanisms: formation of hadrons from gluons, conversion 
of the non-perturbative vacuum energy into particles, or hadronization of a quark gas at 
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constant temperature. Hadronization of a quark gas should approximately conserve the net 
number of particles due to the constraint of entropy conservation. Since hadrons are formed 
of two or more quarks, creation of quark-antiquark pairs should accompany hadronization. 
All three mechanisms for late-stage quark production involve a change in the degrees of 
freedom. Therefore, any signal that pinpoints the time where quarks first appear in a colli- 
sion would provide valuable insight into understanding whether a novel state of matter has 
been formed and persisted for a substantial time. 

The link between balance functions and the time at which quarks are created has a 
simple physical explanation. Charge-anticharge pairs are created at the same location in 
space-time, and are correlated in rapidity due to the strong collective expansion inherent to 
a relativistic heavy ion collision. Pairs created earlier can separate further in rapidity due 
to the higher initial temperature and due to the diffusive interactions with other particles. 
The balance function, which describes the momentum of the accompanying antiparticle, 
quantifies this correlation. The balance function describes the conditional probability that 
a particle in the bin p\ will be accompanied by a particle of opposite charge in the bin p 2 : 

B{pi\P\) = -{p(b,p 2 \a,p 1 ) - p(b,p 2 \b,p 1 ) +p(a,p 2 \b,pi) - p(a,p 2 \a,pi)} , (2) 

where p(b,p2\a,pi) is the conditional probability of observing a particle of type b in bin p 2 
given the existence of a particle of type a in bin p\. The label a might refer to all negative 
kaons with b referring to all positive kaons, or a might refer to all hadrons with a strange 
quark while b refers to all hadrons with an antistrange quark. Balance functions will be 
discussed in greater detail in other articles of these proceedings [|28|, |9| pofl . 

6. Conclusions 

The Q./Q. has been found to be extremely sensitive on the dynamics of hadronization in 
p+p reactions. Within the fragmenting color flux tube models we have predicted that Q./Q. 
ratio is significantly above unity. This is in strong contrast to statistical model results which 
always imply that B/B ratios are bellow or equal to unity in proton+proton reactions. Since 
this observable is accessible by NA49 measurements at the SPS it can provide an excellent 
test to distinguish the statistical model hadronization scenario from that of microscopic 
color-flux tube dynamics - first very preliminary results seem to support the statistical 
hadronization hypothesis. Balance functions yield complimentary information and offer a 
unique, model-independent formalism to probe the time-scale of hadronization in heavy 
nucleus-nucleus collisions. 
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